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La presente invention conceme une methode pour simuler des tests de 
puits sur un modele discret de reservoir fissure. 

5 La methode peut Stre mise en oeuvre par exemple dans le domaine de la 

production p6troliere par des mgenieurs de gisement en vue d'obtenir des 
predictions d'ecoulement fiables. 

n s'agit de simuler un essai de puits dans un milieu poreux permeable 
(reservoir) traverse par un reseau de fissures beaucoup plus conductrices que 

10 la matrice poreuse. Les reservoirs fissures constituent un type extreme de 
reservoirs heterogenes comportant deux milieux contrastes, un milieu matriciel 
contenant la plus grande part de ITiuile en place et presentant une faible 
permeability, et un milieu fissure representant moins de 1% de ITiuile en place 
et hautement conducteur. Le milieu fissure lui-meme peut Stre complexe, avec 

15 differents ensembles de fissures caracterises par leur densite, longueur, 
orientation, inclinaison et ouverture respectives. 

Lors d'un essai de puits, les conditions de debit imposees au puits 
amenent ITiuile contenue dans le reservoir a s'ecouler vers le puits. II s'agit d'un 
ecoulement monophasique (la seule phase mobile est ITiuile) et compressible 
20 (te roche du reservoir et le fluide huile sont compressibles). 

Pour tout volume elementaire du reservoir, la pression de ITiuile 
contenue dans ce volume obeit, si on ne tient pas compte de la gravite, k 
l'equation suivante : 



3P K 
(t>-C T .— = div(-.VP) + Q (1) 

ou: 

<(> :represente le volume poreux, 
Ct, la compressibilite totale (fluide +• roche), 
5 K, la permdabilite de la roche, 

jn, la viscosite du fluide, 

Q, le debit entrant qui est nul partout sauf aux endroits oil le puits est en 
communication avec le reservoir, et 
P, ia pression inconnue. 

10 Pour simuler un essai de puits, dans quelque milieu que ce soit, on doit 

rdsoudre cette Equation dans l'espace et dans le temps. On realise done une 
discretisation du reservoir (maillage) et la resolution du probleme consiste a 
trouver les pressions des mailles au cours du temps, lui-meme discretise en un 
certain nombre de pas de temps. 

15 Dans le cas d'un milieu fissure, si Ton veut tenir compte de toute la 

complexite du reseau de fissures, il est necessaire de realiser un maillage qui 
reprdsente correctement ce reseau. Par ailleurs, dans un tel milieu, les 
permeabilites K sont en general beaucoup plus elevees dans les fissures que 
dans la matrice et les ecoulements sont par consequents rapides dans les 

20 fissures et lents dans la matrice. 

Etat de la technique 



Actuellement on utilise des outils de calcul d T ecoulements 
monophasiques mais ils sont appliques non au milieu geologique reel dans 
toute sa complexite mais a une representation homogeneisee, suivant le 
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modele de reservoir dit de double milieu decrit par exemple par Warren et 
Root, dans « The Behavior of Naturally Fractured Reservoirs », SPE Journal 
Septembre 1963. Tout volume elementaire du reservoir fissure est ainsi 
modelise sous la forme d'un ensemble de blocs parallelepipediques identiques 

5 limites par un systeme orthogonal de fissures uniformes continues orientees 
dans la direction de Tun des trois principaux sens d'ecoulement (modele dit de 
« boite a sucre »). L'ecoulement des fluides a l'dchelle du reservoir s'effectue a 
travers le milieu de fissure seulement et des echanges de fluides surviennent 
localement entire les fissures et les blocs matriciels ne sont pas appliques au 

10 milieu complexe. Cette representation qui ne rend pas compte de la complexite 
du reseau de fractures dans un reservoir, se revele efficace neanmoins mais a 
1' echelon d'une maille de reservoir ay ant typiquement pour dimensions 100m 
x 100m. 

n est bien preferable cependant d'effectuer les calculs d'ecoulement sur 
15 le modele geologique « reel » du miUeu et non sur un modele homogene 
equivalent, car cela presente le double avantage de pennettre : 

- de valider l'image geologique du reservoir construite par le geologue a partir 
de 1'ensemble des informations qu'il a pu recueillir sur la fracturation du 
reservoir (cette validation s'effectuant par comparaison aux donnees d'essais 

20 de puits) ; 

- de tester de maniere fiable la sensibilite du comportement hydraulique du 
milieu aux incertitudes sur l'image geologique du milieu fissure. 

Une methode de modelisation connue consiste a mailler finement le 
reseau de fissures et la matrice en ne faisant aucune approximation concernant 
25 les echanges de fluides entre les deux milieux. Elle est cependant difficile a 




mettre en ceuvre car la geometrie souvent complexe des espaces entre les 
fissures rend difficile leur maillage et qu'en tout etat de cause, le nombre de 
mailles a traiter finalernent est souvent enorme. La complexite augmente 
encore s'il s'agit d'un maillage en 3D. 

5 Des techniques de modelisation de milieux poreux fissures sont decrites 

dans les brevets FR-A-2.757.947 et FR-A-2.757.957 . du demandeur. La 
premiere est relative a la determination de la permeabilite de fissure 
6quivalente d'un reseau de fissures dans un milieu multi-couche souterrain a 
partir d'une representation connue de ce reseau, permettant de relier de 

10 maniere systematique des modeles de caracterisation de reservoir fissure a des 
simulateurs double porosite en vue de realiser une modelisation plus realiste 
d'une structure geologique souterraine fissuree. 

La deuxieme est relative a la modelisation simplifiee d*un milieu 
geologique heterogene poreux (tel qu'un reservoir traverse par un reseau 
15 irregulier de fissures par exemple) sous la forme d'un milieu transpose ou 
equivalent de maniere a ce que le milieu transpose soit equivalent au milieu 
d'origine, relativement a un type determine de fonction de transfert physique 
(connu pour le milieu transpose). 

La methode selon T invention permet de simuler les ecoulements 
20 monophasiques dans un milieu fissure. 

La methode selon T invention a pour objet pour modeliser les flux de 
fluides dans un milieu poreux multi-couches fissure en tenant compte de la 
geometrie reelle du reseau de fissures et des echanges locaux entre la ma trice 
poreuse et les fissures en chaque noeud du reseau et ainsi de permettre une 
25 simulation des interactions entre les variations de la pression et du debit au 



'p. 
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niveau d'un puits traversant le dit milieu. Elle est caracterisee en ce que Ton 
discretise le milieu fissure par un maillage, dans lequel les mailles de fissure 
sont centrees sur des noeuds aux differentes intersections des fissures, chaque 
noeud etant associe a un volume de matrice, et Ton deterrrune les flux entre 
5 chaque maille de fissure et le volume de matrice associe dans un regime de 
flux semi-permanent. 

On delimite par exemple le volume de matrice associe dans chaque 
couche, a chaque maille de fissure, par 1* ensemble des points plus proches du 
noeud correspondant que de noeuds voisins. 

10 Pour determiner le volume de matrice associe a chaque maille de 

fissure, on discretise par exemple chaque couche fissuree en pixels et on 
determine la distance de chaque pixel a la maille de fissure la plus proche. 

La valeur de la transmissivite pour chaque couple maille de fissure - 
bloc matriciel est determinee par exemple en considerant que la pression varie 
15 lineairement en fonction de la distance du point considere a la maille de fissure 
associee au bloc. 

La methode selon F invention simplifie beaucoup le calcul des echanges 
entre la matrice et les fissures. En effet, il n'y a pas de maillage specifique de 
la matrice mais chaque noeud du reseau de fissure est associe a un volume de 
20 matrice, et l'ecoulement entre chaque noeud fissure et son bloc associe est 
calcule en regime semi-permanent. L'innovation concerne le calcul du volume 
de chaque bloc ainsi que du coefficient de proportionalite (appele 
transmissivite) qui relie de debit matrice- fissure a la difference de pression 
matrice- fissure. 
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• L'avantage majeur de la methode selon 1* invention par rapport k 
celles operant sur reseau fissure reel flnement maille, tient, comme on le verra 
mieux dans la suite de la description, a la reduction considerable du nombre de 
mailles centrees sur des noeuds de fissure utilisees pour la resolution de 
5 1' equation (1). Le gain en vitesse de simulation par rapport aux methodes 
anterieures k maillage fin, on peut le verifier, est supdrieur au simple rapport 
arithmetique du nombre de mailles utilise dans les deux cas. 

La methode selon 1' invention D'autres caracteristiques et avantages de 
la methode selon rinvention, apparaitront k la lecture de la description ci-apres 
10 d'un exemple non limitatif de realisation, en se referant aux dessins annexes 
ou: 

- la Fig.l montre un exemple de maille de fissure en 2D ; 

- la Fig. 2 montre un exemple de bloc matriciel obtenu en 2D ; 

- la Fig. 3 montre une courbe de variation de debit imposee dans un essai de 
15 puits ; et 

- la Fig. 4 montre un exemple de reseau de fissures pour lequel la methode 
selon rinvention conduit a une reduction radicaie du nombre de mailles a 
traiter. 

DESCRIPTION DE LA METHODE 

20 1) Discretisation du reseau de fissures 

La methode comporte un maillage du reservoir fissure multi-couches tel 
qu'il est modelise en utilisant par exemple la technique de modelisation de 
reservoirs fissures decrite notamment dans le brevet FR 2.757.947 precite, en 
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admettant que les fissures sont sensiblement perpendiculaires aux limites des 
couches. 

Pour former ce maillage, on repere toutes les intersections de fissures 
ou noeuds (Fig.l) et on constitue des mailles de fissure en associant a chaque 

5 noeud un bloc matriciel unique. Pour les besoins d'une modelisation en 3D, on 
est conduit a ajouter des noeuds supplementaires dans les couches voisines 
afin de tenir compte des fissures qui en traversent plusieurs. Les noeuds de 
calcul ainsi definis constituent le centre des mailles de fissure. On fixe comme 
limites aux mailles, les extremites des fissures d'une part et le milieu des 

10 segments reliant les noeuds de calcul d'autre part. Compte-tenu de ces limites 
imposees aux mailles, on peut calculer leur volume $ (cf. relation 1). Le calcul 
des liaisons entre mailles de fissure (transmissivites) utilisees pour le calcul 
des flux inter-mailles peut Stre effectue selon la methode decrite dans le brevet 
2.757.947 precite. : 

15 2) Discretisation du milieu matriciel 

La discretisation du milieu matriciel consiste k affecter un volume de 
roche h. chacune des mailles de fissure. Lors de la simulation dynamique, les 
echanges entre la matrice et les fissures seront calcules en mode semi- 
permanent (flux d'echange proportiormel a la difference de pression) entre 
20 chaque maille fissure et son unique bloc matriciel associe. 

Pour le calcul de ces volumes de matrice, on traite le probleme couche 
par couche, ce qui revient a negliger les ecoulements verticaux par rapport aux 
ecoulements horizontaux dans la matrice. 

Pour une couche donnee, on definit les blocs de la maniere suivante : 
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- les hauteurs de blocs sont egales et fixees a la hauteur de la couche ; 

- dans le plan de la couche, la surface du bloc associ6 a une maille fissure est 
l'ensemble des points plus proches de cette maille fissure que d'une autre. 

Physiquement, cette definition suppose que les limites a flux nul dans la 
matrice sont les lignes d'equidistance entre les mailles de fissure. Cette 
approximation est valable si les fissures voisines sont a des pressions tres 
proches, ce qui est vrai lorsqu'on est en presence de fissures fortement 
conduc trices par rapport a la matrice. 

a) Geometrie des blocs 

Pour detenniner la geometrie des blocs dans une couche donnee, il 
s'agit done de resoudre un probleme bi-dimensionnel qui consiste a trouver, 
pour chaque maille fissure les points les plus proches de cette maille que d'une 
autre. Ce probleme est resolu en utilisant par exemple mais de preference la 
methode geometrique exposee dans 1' autre brevet FR 2.757.957 precite, qui 
pennet, en discretisant la couche fracturee en un ensemble de pixels et en 
appliquant un algorithme de traitement d'image, de determiner la distance de 
chaque pixel a la fracture la plus proche. Dans la phase ^initialisation de cet 
algorithme, on affect la valeur 0 (distance nulle) aux pixels appartenant a une 
fracture et une grande valeur aux autres. Si dans cette phase, on donne, de 
plus, pour chaque pixel de fissure, le numero de maille de fissure auquel il 
app artient, le meme algorithme pennet finalement de detenniner, pour chaque 
pixel : 

• la distance de ce pixel a la maille fissure la plus proche ; 

• le numero de la maille fissure la plus proche. 
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L'ensemble des pixels ayant le meme numero de maille de fissure 
constitue la surface du bloc matriciel associe a cette maille. En multipliant 
cette surface par la hauteur de la couche, on obtient le volume du bloc associe 
a la maille. La Fig.2 montre un exemple de bloc matriciel ainsi obtenu en 2D : 

5 Par construction, la somme des surfaces des blocs est egale a la surface 

totale de la couche, ce qui garantit qu'on ne rajoute ni n'enleve du volume a la 
couche. 

Pour chaque bloc matriciel, ralgorithme de traitement d'image donne 
aussi la distance de chaque pixel du bloc a la maille fissure associee. On se 
10 sert de cette information pour calculer la transmissivite entre la maille fissure 
et le bloc matriciel. 

b) Transmissivite matrice-fissure 

Dans llrypothese d'un regime semi-permanent ou Ton considere que le 
flux d'echange est proportionnel a la difference des pressions de Tun et l'autre, 
15 on a la relation suivante : 

Fmf = I J L -(Pm-Pf) (1) 

ou: 

Fmf est le flux matrice-fissure, 

Tmf, la transmissivite matrice-fissure, 
20 la viscosite, 

P m , la pression du bloc matriciel, et 

P f , la pression de la maille fissure associee. 
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Dans lliypothese du regime semi-permanent, la valeur de la 
transmissivite Tmf est constante au cours de la simulation. Cette valeur est ici 
calculee pour chacun des couples maille de fissure - bloc matriciel. 

Pour ce calcul, on fait l'hypothese que dans le bloc matriciel, la pression 
5 varie lineairement en fonction de la distance du point considere a la maille 
fissure associee au bloc. On definit la transmissivite comme suit : 



1 est la longueur des fissures dans la maille fissure, 
10 H, la hauteur de la couche (et du bloc matriciel), 
K, la permeabilite de la matrice 

D, la distance des fissures a laquelle la pression est la pression moyenne du 



On connait 1, H et K (le facteur 2 vient du fait que la fissure a deux 
15 faces). Le calcul de D est fait a partir des informations donnees par 
ralgorithme de traitement d'image sur les distances des pixels aux fractures les 



T mf = 



21-H-K 
D 



(2) 



ou: 



bloc 



plus proches. En effet, suivant 1'hypothese de variation lineaire de la pression 
en fonction de la distance a la fissure, on a : 




(3) 



20 



ou S est, en 2D, la surface du bloc matriciel et d(s) est la fonction donnant la 
distance aux fissures des points du bloc matriciel a la maille fissure associee. 
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En terme de pixels, on a done simplement : 

D = ^-Xdn (4) 
iN N 

ou N est le nombre de pixels du bloc matriciel et d n la distance du pixel n a la 
maille fissure. 

5 3) Representation des puits 

Un puits est represents par sa geometrie d'une part et les contraintes de 
d6bit qui lui sont imposdes au cours du temps d'autre part. 

a) Geometrie 

Un puits est constitue par un ensemble de segments connectes qui 
10 intersectent les fractures du reseau. La representation geometrique d'un puits 
est done une ligne brisee en 3D. Certains segments de cette ligne brisee sont 
en communication avec le reservoir. Les points de communication entre le 
puits et le reseau sont les points d'intersection de ces segments communicants 
avec les fractures. 

15 b) Contraintes de debit 

Les contraintes de debit sont imposees au point ou le puits entre dans le 
milieu fracture. Elles sont entrees par l'utilisateur sous la forme d'une fonction 
en plateau donnant le debit en fonction du temps. Les temps de changement de 
debit du puits indiquent les differentes periodes a simuler. 

20 Pour un essai de puits classique (Fig. 3), on impose un debit pendant une 

premiere periode au terme de laquelle on ferme le puits (debit nul) et Ton 
observe la remontee de pression dans le puits. 
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4) Simulation dynamique 

Lors de la simulation dynamique, la periode de temps simulee est 
decomposed en pas de temps dont la duree est comprise pour un essai de puits 
entre 1 s (juste apres rouverture ou la fermeture du puits) et quelques heures. 

Le passage de l'instant n (auquel on connait toutes ies valeurs de 
pressions) et l'instant n+1, est realise en resolvant dans toutes les mailles et 
blocs du reservoir, l'equation suivante qui correspond a la relation 1 
discretisee : 

p n+l_ p n T:- 

t>i-CT-^ Tr 1 ^ + i:f-(pr 1 -pr 1 ) = Qi e> 

I I J 

ou : 

i, est le numero de la maille ou du bloc, 
j, le numero de la maille ou du bloc voisin de i, 
P n , la pression a l'instant n, 
t n , le temps a rinstant n, 
Qi, , le debit entrant dans i, et 
Ty , la transmissivite de liaison entre i et j. 

A part les pressions, les termes de cette equation sont connus. Les 
volumes poreux des mailles de fissures et des blocs matriciels (<&) sont connus 
par le maillage, les transmissivites fissure-fissure et matrice- fissure (Tij) sont 
calculees comme decrit plus haut et les debits (QO sont nuls partout sauf aux 
connexions puits-reservoir ou ils sont imposes. 
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a) Elimination des inconnues relatives a la matrice 

L'equation (6) peut etre resolue par l'inversion d'un systeme lineaire ou 
les inconnues sont les pressions des mailles de fissure et les pressions des bloc 
matriciels a l'instant n+1. 

Neanmoins, comme chaque maille de matrice n'est reliee qu'a un seul 
bloc matriciel et que les blocs ne sont pas relies entre eux, il est possible 
d'elirniner les inconnues des blocs du systeme dont la taille est alors divise par 
2. 

En effet, pour un bloc matriciel, l'equation (6) s'exprime comme suit : 

<t>m*C T m- t n+l_ t n + ^ " (P ^ " Pf } " ° W 

les indices met f signifiant respectivement matrice et fissure, et Cxm 
etant la compressibilite totale de la matrice. 

On a C Xm =c m + Sw-c w +(l-Sw)-c 0 , ou Sw est la saturation 
residuelle en eau dans la matrice, c m la compressibilite de la matrice, c w la 
compressibilite de Teau et c Q la compressibilite de lliuile, et (|> m est le volume 
poreux du bloc matriciel, c'est-a-dire le volume du bloc multiplie par la 
porosite de la matrice. 

On en deduit : 
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avec : 



A 4>m' C Tm 
A m - t n+l _ t n 



En reportant cette expression dans 1'equation concernant la maille de 
fissure, on obtient : 



m m 



ou i est le numero de la maille de fissure consideree, j, les numeros des mailles 
fissure connectees a i et m, le numero du bloc associe a la maille de fissure i. 

De plus, on a Ai = ^\F'' TS ou C T f est la compressibility totale de 
fissure : C Tf = c f + c 0 . 

La reduction d'un facteur 2 du nombre d'inconnues est un avantage de 
plus de la presente methode. Le temps de resolution d'un systeme lineaire 
evolue . approximativement comme le carre du nombre d'inconnues, par 
consequent, Felimination des inconnues matrice entrame un gain important 
dans la vitesse de simulation d'un essai de puits sur un reseau fissure. 

b) Algorithme de resolution d'un pas de temps 

L'algorithme general de resolution d'un pas de temps consiste a 
consumer un systeme lineaire correspondant a Tequation (9) ecrite pour toute 
maille fissure i, a resoudre ce systeme lineaire et a mettre a jour les inconnues. 
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Les etapes de cet aigorithme sont les suivantes : 

- Constitution du systeme lineaire avec calcul des termes de l'equation (9) 

Calcul de la contribution du milieu fissure 
tenne d'accumulation (AO 
5 liaison entre mailles fissure (Tq/^) 

conditions limites puits (QO 

Calcul de la contribution du milieu matrice 
terme d'accumulation (Am) 
tenne de liaison matrice fissure (UJ 

10 - Resolution du systeme lineaire (pour trouver les valeurs de pression de 
fissure & Tinstant n+1) par l'algorithme iteratif du gradient conjugue (CGS) 

- Mise k jour avec le calcul des pressions des blocs matriciels en utilisant 
l'equation (8), et gestion du pas de temps. 

c) Gestion du pas de temps 

15 Dans d'une simulation d'essai de puits, les pas de temps sont en general 

tres courts apres rouverture ou la fenneture du puits, parce que c'est a ces 
moments que les variations de pressions sont fortes. Plus tard, les pas de temps 
peuvent etre plus longs, lorsque les variations de pression sont plus faibles. 

La gestion du pas de temps consiste a relier la duree du pas de temps 
20 n+1 k la variation de pression maximum observee lors du pas de temps n. Pour 
ce faire, 1'utilisateur doit fournir les donnees suivantes : 

Pas de temps initial : dtini 
Pas de temps minimal : dt^ 



16 



Pas de temps maximum : dtmax 
Borne inferieure de variation de pression : dpnun 
Borne superieure de variation de pression : dpmax 
Facteur d'augmentation du pas de temps : rdt+ (>1) 
Facteur de reduction du pas de temps : rdt- (<1) 

A partir de ces valeurs, la gestion du pas de temps est r6alise*e comme 

suit : 

A rinstant t n , on calcule le majorant des variations de pressions dans les 
mailles fissure et les blocs matriciels pendant le pas de temps n (c'est a dire 
entre les instants f~ l et t n ). Suivant la valeur de cette variation de pression dp, 
le pas de temps est soit : 

• augmente d'un rapport rdt+ si dp < dpmin, 

• diminue d'un rapport rdt- si dp > dp 

max? 

• inchange dans les autres cas. 

Apres chaque changement de condition de debit au puits, le pas de 
temps est reinitialise a la valeur dt^. De plus, une optimisation est realisee en 
fm de periode de simulation : si tf est le temps a simuler pour atteindre la fin 
de periode et dt le pas de temps prevu, le pas de temps choisi est min(t f /2,dt). 

6) Donnees d'entree de la simulation 

a) Reseau de fractures 

La geometrie du reseau de fractures ainsi que les attributs des fractures 
(conductivite, ouverture) sont donnes sou la forme d'un fichier comme dans la 
methode decrite dans le brevet FR 2.757.947 precite. 
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b) Petrophysique 

Les donnees p6trophysiques a fournir sont les suivantes (on indique 
runite entre parenthese) : 

• compressibilite du reseau de fissures (c f en bar" 1 ) 
© compressibilite de la matrice (c m en bar' 1 ) 

• permeabilite de matrice (K en mD) 

• porosite de matrice (4> m sans dimension) 

Les donnees en question sont considerees comme homogenes sur le 
milieu fissure considere. 

c) Fluide 

Les donnees concemant les fluides contenues dans le milieu fissure 
concernent l'huile (mobile) et Teau (immobile) qui est toujours presente en 
quantite residuelle dans la matrice : 

• saturation residuelle en eau dans la matrice (Sw sans dimension) 

• compressibilite de l'eau (c w en bar" 1 ) 

• compressibilite de l'huile (c 0 en bar* 1 ) 

• viscosite de l'huile (ja en cpo) 

d) Puits 

Les donnees concemant le puits sont : 

• sa geometrie sous la forme d'une serie de segments connectes 



18 




• les debits imposes sous la forme (Tune courbe dormant le debit 
impose en fonction du temps. 

e) Gestion du temps 

Les valeurs a fournir sont les valeurs deja definies dans le chapitre sur 
5 la gestion du pas de temps : dt^, dW, dt^, dpmin, dpmax, rdt+ et rdt 

Exemple comparatif de maillage 

L'exemple comparatif suivant illustre bien la simplification du maillage 
et le gain correlatif en temps de calcul des flux, que Ton peut obtenir par la 
mise en ceuvre de la m6thode. A titre de comparaison, on a observe que la 
10 zone fissuree representee a la Fig.4, que Ton aurait maille normalement avec 
80.000 mailles avec les techniques de maillage classiques, peut etre maille ici 
avec seulement quelque 4000 mailles. 
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REVINDICATIONS 



1) Methode pour modeliser les flux de fluides dans un milieu poreux 
multi-couches fissure dans le but de simuler les interactions entre les 
variations de la pression et du debit au niveau d'un puits traversant le dit 
milieu, caracterisee en ce que Ton discretise le milieu fissure par un maillage, 
dans lequel les mailles de fissure sont centrees sur des noeuds aux differentes 
intersections des fissures, chaque noeud etant associe a un volume de matrice, 
et Ton determine les flux entre chaque maille de fissure et le volume de 
matrice associe dans un regime de flux semi-permanent. 

* 2) Methode selon la revendication 1, caracterisee en ce que le volume 
de matrice associe a chaque maille de fissure est delimite dans chaque couche 
par F ensemble des points plus proches du noeud correspondant que de noeuds 
voisins. 

3) Methode selon la revendication 2, caracterisee en ce que Fon 
discretise chaque couche fissuree en pixels et on delimite le volume de matrice 
associe h chaque maille de fissure par une detennination de la distance de 
chaque pixel a la maille de fissure la plus proche. 

4) Methode selon Fune des revendications precedentes, caracterisee en 
ce que Fon determine la valeur de la transmissivite pour chaque couple maille 
de fissure - bloc matriciel en considerant que la pression varie lineairement en 
fonction de la distance du point considere a la maille de fissure associee au 
bloc. 
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